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Photo-induced forces can be used to manipulate and cool the mechanical motion of oscil¬ 
lators. When the oscillator is used as a force sensor, such as in atomic force microscopy, 
active feedback is an enticing route to enhancing measurement performance. Here, we show 
broadband multimode cooling of —23 dB down to a temperature of 8 ± 1 K in the stationary 
regime. Through the use of periodic quiescence feedback cooling, we show improved signal- 
to-noise ratios for the measurement of transient signals. We compare the performance of real 
feedback to numerical post-processing of data and show that both methods produce similar 
improvements to the signal-to-noise ratio of force measurements. We achieved a room tem¬ 
perature force measurement sensitivity of < 2 x 10 16 N with integration time of less than 
0.1 ms. The high precision and fast force microscopy results presented will potentially benefit 
applications in biosensing, molecular metrology, subsurface imaging and accelerometry. 

m 

Introduction 

To attain high-resolution mass, position and force sensing, cantilevers with low effective masses 
and high mechanical quality factors are required. Nano-scale resonators are excellent candi- 
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dates for mass and force sensing 00. biomaterial sensing by means of Kelvin probe force mi¬ 
croscopy (4j, and single-spin and charge detection [5j. Miniaturized oscillators, however, are more 
susceptible to thermal noise and their use can be challenging as it entails low-noise high-frequency 
electronics to monitor and control the vibrational modes. Passive cooling is not always an option, 
particularly in biomaterial sensing applications where samples are dependent on specific environ¬ 
mental conditions. Active feedback cooling can instead be used to reduce the Brownian motion of 
the oscillators. 


To date, effective optical cooling of micro/nano-mechanical resonators has been achieved 
via cavity cooling [|6]-|8j and active optical feedback cooling (9j through radiation pressure [(T0|-[T3 1 


and photothermal |14l 151 forces. A feedback system will suppress all motion of the oscillator and 
has no way to differentiate between motion due to thermal noise and a signal that one wishes to 
detect. The signal-to-noise ratio (SNR) of the sensor under conditions of steady-state feedback 


therefore may not be straightforwardly enhanced [16]. Not all is lost, however, as the increased 


effective damping rate due to feedback can give the oscillator a larger linewidth, which means that 
the oscillator will respond more quickly to external forces. As a result, less integration time is 
required to detect steady-state signals with the adoption of feedback cooling [|T7|. 


Feedback cooling requires constant measurement of the system in order to generate a control 
signal to counteract thermalisation. If this measurement data is not used for feedback, one may 
ask if the data can still be used to improve our knowledge of the system dynamics and the mea¬ 
surement outcome. Provided that we are dealing with linear systems, this measurement data can 
indeed be fed into numerical estimation algorithms such as Kalman or Wiener filters. These tech¬ 
niques allow tracking and prediction of the motion, which can yield advantages similar to those 
of feedback cooling. In fact, it is possible to use numerical algorithms to simulate the effect of 


physical feedback cooling [18]. These estimation techniques require accurate knowledge of the 


system parameters to perform efficiently. In a nano-scale opto-mechanical system, the phase and 


amplitude noise of the laser can alter the system parameters [ 191 and real-time tracking may be 
required for optimum estimation. 
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Outside the confines of the steady-state regime, feedback and estimation can be used to im¬ 


prove the SNR [16 18]. To date, active cooling of nano-mechanical systems in the transient regime 
has not been observed and its implications for impulsive force sensing have not been studied. In 
the presence of an impulsive force, feedback cooling can be turned on prior to measurement and 
then turned off during the measurement. This periodic quiescence feedback works in such a way 
that the effect of an impulsive external force on the resonator is measured before the system is fully 
thermalized. In this case, the SNR of a pulsed signal can be improved by synchronizing the feed¬ 
back in cycle with the signal. As with the steady-state conditions, it has also been predicted that 
an equivalent enhancement can be achieved by using estimation techniques provided the system 


dynamics are well known (see Supplementary Information of Ref. 18]). 


In this letter we investigate the dynamics of gold-coated Ag 2 Ga nanowire^J Using homo¬ 
dyne detection we can measure the motion of the nanowire and then actively control the vibration 
of its modes using bolometric forces. Our setup is cavity-free, allowing cooling to be accom¬ 


plished over a wider bandwidth where simultaneous cooling of multiple vibrational modes [20] 
is observed. We explore the implementation of non-stationary active cooling on nano-mechanical 
oscillators by means of optical forces and observe SNR enhancement of impulsive force measure¬ 
ments. We then compare the results with off-line estimation techniques and show enhancement by 
a factor of about 5 using periodic feedback cooling as well as estimation methods. Using both phys¬ 
ical control and estimation techniques we achieved a force sensitivity of better than 2 x 10 -16 N. 


Results 

Detection and Feedback. The experimental setup is presented in Fig.[l}a), with SEM images of 
a nanowire shown in Fig. jT;b) (for more details on its properties and characterization refer to Sup¬ 
plementary Note 1). Optomechanical cooling often relies on optical cavity resonators to enhance 
the coupling to the mechanical modes. In that case the operational bandwidth is limited by the 
cavity dynamics and allows only a single mechanical frequency to be addressed via optical inter- 
1 http://nauganeedles.com 
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action. In our experiment, we use a microscope objective to maximize the single-pass interaction 
of the nanowire with the laser light. The absence of a cavity allows simultaneous access to all of 
the nanowire’s mechanical resonances. Light scattered from the nanowire is used in a homodyne 
measurement system j9[. The motion of the nanowire is therefore referenced to the phase of the 
detection laser. The signal-to-noise ratio of this measurement is limited by the amount of scattering 
from the nanowire. The power used in the measurement is low enough to prevent the observation 
of any back-action due to the detection laser. 


The actuation laser has a different wavelength to the detection laser to avoid any chance of 
interference between these subsystems. The dominant force on the nanowire is bolometric, i.e. it 


results from differential thermal expansion of a bimorph structure [SI -23]. In our case, it is the 
composite Au/Ag 2 Ga structure of our nanowires that leads to strong thermal effects. The speed of 
the force depends on the thermal time constant of the object being heated, so it is naturally slow 
for macroscopic oscillators. For nanoscale objects, however, it can be both fast and substantial 


in magnitude [24 25). The response time of the nanowire to laser driving is found to be 15 ps 
at ambient pressure. Unlike radiation pressure, the direction of the force depends only on the 
structure of the object and is independent of the direction of the incident beam that is used to drive 
the nanowire. Experimental evidence of these aspects is offered in Supplementary Notes 2 and 3 
and Supplementary Figures 1, 2, 3 and 4. Bolometrically actuated feedback cooling is somewhat 
counterintuitive, as the amplitude of the oscillations is subdued by use of a driving force of thermal 
origins - effectively cooling by heating. Even so, recent theoretical work shows that photothermal 


effects can assist cooling towards the quantum ground state [26-28 [ and photothermal back action 


in a cavity has also been used to cool a semiconductor membrane [ 151 and a graphene sheet 14]. 


Fig. [2} a) shows the displacement spectrum of the nanowire for increasing power of the feed¬ 
back beam in vacuum conditions. The phase of the feedback is adjusted using an adjustable band¬ 
pass analog filter. The feedback phase can be tuned such that simultaneous cooling of multiple 
modes of the nanowire can be obtained, as seen in Fig. |2}b) where we plot the uncooled (red) and 
cooled (blue) vibrational modes of the nanowire up to 2 MHz for a particular phase and gain of 
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the feedback loop. Applying a more advanced RF phase shifter, it is in fact feasible to efficiently 
cool vibrational modes over a broader frequency range. A comparison of the best effective tem¬ 
peratures reached in ambient and vacuum conditions is plotted as a function of feedback strength 
in Fig. gc); we have achieved cooling greater by an order of magnitude and for reduced driving 
power in vacuum compared to ambient pressure due to the improved mechanical quality factor. 
Fig. gd) shows the relation between phase and cooling for the first two oscillation modes of a 
nanowire in ambient pressure conditions; the difference in efficiency is attributed to dissimilarities 
in geometrical qualities along the directions of oscillation. 

SNR Enhancement. To measure the SNR in the transient regime, a 0.1 ms long signal modulated 
at the mechanical frequency is sent to the nanowire right after the feedback is turned off (see 
Fig. (TJa) and (c)). At the moment the signal arrives, with the feedback just switched off, the 
nanowire will also begin to thermalise at the rate of its mechanical dissipation, 7 m . Provided the 
signal time (r s ig ) is much shorter then the re-thermalisation time (7“ 1 ), the measurement of the 
signal will not be significantly affected by thermal noise. In vacuum conditions the mechanical 
dissipation rate of the nanowire is lower than 1 kHz, allowing integration times up to 1 ms long. 
We integrate the energy of the signal after introducing the impulsive force and compare it with the 
thermal noise to estimate the SNR in the presence and absence of periodic cooling. 

As we outlined earlier, linear feedback can also be simulated off-line using estimation tech¬ 
niques. We have performed both Kalman filtering and virtual feedback cooling (see Methods) on 
data taken in the absence of feedback. The digital Kalman filter is constructed based on prediction 
and update stages. At the prediction stage, the model of the system is used to estimate the evolution 
of the system for a short time into the future. In the update stage, the actual measurement results 
and known measurement and process noise vectors are used to refine the estimated evolution. The 
update stage of the filter is switched off when the impulsive force arrives, and the sizes of signal 
and noise are calculated as the phase-space distance between the measured and estimated trajec¬ 
tories in the presence and absence of the impulsive force. The virtual cooling method is used to 
directly simulate ideal periodic feedback cooling. A comparison of the estimation techniques with 
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physical cooling in the transient regime is shown in Fig. [3] Physical feedback cooling is shown to 


be as effective as the virtual cooling, which indicates near optimal actuation of the nanowire. For 
integration times longer than 0.4 ms, feedback cooling shows slightly higher improvement than 
virtual cooling. This is likely due to perturbation of the system parameters by laser noise within 
the time scale of mechanical decay time. The Kalman filter is shown to outperform both physical 
and virtual cooling. In all cases, the SNR degrades with a rate corresponding to the mechanical 
decay time as the thermal noise contribution increases. The SNR enhancement factor r] s N r, defined 
as the ratio between the estimated SNR and that of the raw data where no cooling or filtering is 
used, is plotted in Fig. [3jb). An enhancement factor of about 5 has been observed for both the 
virtual cooling and Kalman filter methods. 

We note here that while estimation methods can provide similar or even higher SNR enhance¬ 
ment factors than feedback cooling for short integration times, they require precise knowledge of 
system dynamics and parameters, and the production of the desired outcome can become compu¬ 
tationally expensive. As an example we noticed that changing the mechanical frequency in the 
model by even 0.1% could significantly change the SNR. Such a change in frequency can occur 
easily through a change in the bulk temperature of the oscillator. 

To calculate the minimum force that can be measured with our setup, we consider the re¬ 
sponse x(t) to a monochromatic force with magnitude F 0 , duration r s j g and frequency ujq/2it that 
is x(t) = 4= J^ sie dt'x(t — t') smfccoF), where y(f) is the mechanical susceptibility. The mean 
squared displacement (for integration time r) is ( x 2 ) = \ \f^' g dt'x{t — t') sin^of 7 ) | 2 , so, 

the force required to obtain a given mean squared displacement is 



( 1 ) 


This suggests a sensitivity as low as 2 x 10 -16 N after integration time of less than 0.1 ms for data 
filtered by the Kalman method. After 0.1 ms the signal is switched off, so the relative contribution 
from the thermal noise increases and the sensitivity degrades. The high resolution and fast force 
sensing are important factors for bio-sensing where conformational changes typically happen on 
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micro- to millisecond time scales [29]. 


Discussion 


In summary, we have demonstrated simultaneous active cooling of fundamental and higher-order 
modes of a bimetallic nanowire by means of optically induced thermal forces, obtaining temper¬ 
atures as low as 8 ± 1 K. In the transient regime, we investigated the effect of feedback cooling 
on the detection of an impulsive signal to demonstrate enhancement in measurement sensitivity. 
Although virtual feedback and estimation theories applied to off-line data provide similar improve¬ 
ments in sensitivity, they require that the system parameters are well known and constant, which 
may be a limitation in practical sensing applications. These technologies are increasingly impor¬ 
tant for sensing applications ranging from monitoring intermolecular interactions p9| to subsur¬ 
face imaging (!2), and generally provide an alternative to passive cooling for systems that can not 
be refrigerated. Both active feedback and numerical post-processing methods can be used to obtain 
an improvement in sensitivity over short time windows, especially advantageous when the system 


dynamics are changing rapidly and long integration times are not accessible [30]. 


Methods 


Feedback cooling. The power spectrum of the measured displacement of the nanowire driven by 


thermal forces in the presence of feedback is [12] 

2 7m To 


S T f r. () 


m 2 - UJ 2 ) 2 + (1 + g) 2 7^ce 2 ’ 


( 2 ) 


where g is the effective gain that depends on the electronic gain, the mechanical quality fac¬ 
tor, the thermal relaxation time and the induced mechanical rigidity resulting from photothermal 


forces [221, uj m /is the mechanical frequency, typically around 300 kHz, 7 m /27r is the mechani¬ 


cal linewidth, ranging from around 10 kHz in air down to 0.8 kHz in vacuum, m is the nanowire’s 
effective mass, around 5 fg, k B is the Boltzmann constant and T 0 is the initial temperature. 

The vibration of the nanowire introduces phase modulation on the reflected light. We send the DC 
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homodyne signal to a piezo-actuated mirror to lock the phase of the local oscillator beam to that of 
the signal reflected from the nanowire. Above the locking bandwidth, which was around 10 kHz, 
the error signal then becomes a readout of the nanowire motion. The effective temperature and 
mechanical damping rates when feedback is sent to the nanowire are given by T efr = T 0 /(1 + g) 
and 7 eff = (1 + g) 7 m , respectively. The displacement spectrum has to be calibrated to compensate 
for the conversion and operational efficiency of the measurement devices. This is done by the 
integration of the area beneath the spectrum at room temperature and comparing it with the result 
expected from the equipartition theorem. 


For higher gains a correction term [3JJ needs to be applied to the power spectrum (see Supple¬ 
mentary Note 4) and the effective temperature inferred. The latter is then characterized by a lower 
bound, T min = ^/ '''"(y m7 -,S dct , which depends on the spectral density of the detection noise S det . 
The lower the measurement noise and the mechanical linewidth are, the lower the minimum tem¬ 
perature achievable will be. Once this limit is reached, increasing the gain further results in noise 


squashing [31,32] (see Supplementary Note 4 and Supplementary Figure 5). 


Data processing. To measure the SNR, four sets of homodyne signals are recorded at a 
rate of 25 MS s” 1 , in the presence and absence of feedback and in the presence and absence of 
the impulsive force. We then apply a spectral filter to restrict the signal to a 40 kHz bandwidth 
around the mechanical frequency. To calculate the SNR of, for example, feedback cooling data, 
we integrate the energy of the homodyne signal corresponding to the data taken in the presence of 
both periodic feedback and the impulsive force, and divide it by the average integral of the energy 
corresponding to data taken in the presence of periodic feedback and absence of the impulsive 
force. Analogous methods are used to measure SNR in the absence of feedback cooling, and 
when estimation techniques are used in place of physical cooling. All results are averaged over 
150 traces, with the error calculated using standard error (for Fig. [3ja)) or standard deviation (for 

Fig. UN- 


Estimation. We use two different estimation strategies based on the discrete extended Kalman 
filter (EKF) [[33]] and the virtual cooling method [18] to enhance the SNR (refer to Supplementary 
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Note 5 for more details). To identify the state-transition matrix of the EKF, we need to know the 
natural frequency, damping rate, initial amplitude, initial velocity, time interval, process and mea¬ 
surement noise vectors and initial covariance estimates. In order to propagate the state and estimate 
the evolution of the system, we use the Runge-Kutta (or RK4) method to predict the state of the 
system in the next step and update the system’s quadratures according to the relative uncertainties 
of measured and predicted values. We switch off the updating step of the EKF after 1 ms, at the 
time when the impulsive force arrives, so that the EKF predicts the subsequent behaviour of the 
system. This prediction is expected to be reliable while the system is not significantly affected by 
stochastic thermal noise, and deviations from it indicate the presence of the impulsive signal. The 
SNR is calculated in the same way as for feedback cooling, but instead of integrating the energy 
of the oscillator we integrate the phase-space distance between the observed and predicted trajec¬ 
tories. 


The virtual cooling is applied as suggested in the Supplementary Information of Ref. [ 18] in order 
to simulate the periodic feedback cooling scheme. Knowing the mechanical susceptibility of the 
system we perform temporal discretization of the resulting Fredholm equation into a 1000 x 1000 
matrix equation. Solving this yields the simulated measurement record, from which the SNR is 
calculated (refer to Supplementary Note 5 for more details). 
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DETECTION POINT 


(C) 



Figure 1: Schematic experimental setup and experiment cycle, (a) Experimental setup. The 
position of the nanowire is controlled by a nanopositioning stage. One laser (1064 nm) acts as a 
detection beam. Reflected and transmitted light from the nanowire are detected using a balanced 
homodyne and camera, respectively. Light is focused onto, and collected from, the nanowire 
using a pair of microscope objectives (40x, NA = 0.65). A second laser (780 nm) is intensity 
modulated using an acousto-optic modulator, enabling feedback actuation based on the homodyne 
signal. During transient regime measurements, a third laser (795 nm) modulated to provide the 
impulsive signal is used, (b) Mounting arrangement and SEM close-ups of a cylindrical Ag 2 Ga 
nanowire about 60 pm long and 150 nm thick. The detection beam interacts with the free end of the 
nanowire while the feedback or impulsive signal driving beams are typically positioned ~10pm 
from the tip to enhance the thermal bending, (c) Schematic of the experiment duty cycle during 
transient cooling and impulsive force measurement. 
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Figure 2: Experimental results of multimode photothermal cooling of a nanowire, (a) Dis¬ 
placement spectrum of the two orthogonal fundamental modes of a nanowire (R « 100 nm , 
L w 40 pm) for different powers of the feedback beam at vacuum condition (HP 4 Pa), (b) Simul¬ 
taneous cooling of four distinguished resonances of the nanowire in vacuum with phase and gain 
optimized for more efficient cooling on one mode per pair. The displacement spectrum data in the 
absence of feedback is in red and the trace with feedback cooling is in blue. Insets show a close-up 
of the vibrational modes, (c) Effective temperature of the first vibrational mode as a function of 
feedback laser power at ambient (red) and vacuum (blue) conditions. The solid line represents the 
theory assuming a linear relation between laser power and feedback gain g. (d) Effective temper¬ 
ature of the first (dark green) and second (light green) vibrational modes as a function of the phase 
of the feedback signal. Phases of —90° and 90° correspond to maximum oscillation amplification 
and maximum cooling, respectively. The solid lines represent the theory predictions. The error 
bars in (c) and (d) are estimated based on error propagation in the Lorentzian fit of the amplitude 
noise. 
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Figure 3: Comparison of force measurement results of feedback cooling and estimation meth¬ 
ods. (a) Signal-to-noise ratio (SNR) measurement of energy of nanowire driven with impulsive 
photo-thermal force. The black and blue traces show the SNR measurement without and with 
feedback cooling. The red trace is the SNR of non-cooled data after post-processing using Kalman 
filter. The green trace represents the results of virtual cooling on data. The shaded regions represent 
the standard error, (b) The SNR enhancement factor, r/ SNR , defined as the ratio of the estimated 
SNR divided by the SNR of the raw data. The dashed lines represent the standard deviation in 
enhancement estimation, (c) Force resolution as a function of integration time for filtered, raw and 
feedback cooling data. 
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Supplementary Notes 


Supplementary Note 1 - Characterization. The oscillators used in our system are gold-coated 
Ag 2 Ga nanowires (provided by NaugaNeedles LLC [Tj]). They range in size between 20 and 60 pm 
in length and 50 to 200 nm in diameter, and are coated with about 40 nm of gold to improve scat¬ 
tering efficiency, as optical detection of its vibrational modes relies entirely on light scattered 
from its surface. The scattering efficiency (Q sca ) of a long cylinder for light normally incident 
on the nanowire indicates the amount of light scattered from the sub-wavelength object predicted 
by Mie scattering theory [2J. This quantity is plotted as a function of nanowire radius in Suppl. 
Fig. |TJa), with the total scattering efficiency in the XY-plane around the nanowire placed at the 
origin, when laser light is incident from right to left, shown in the inset. It can be used to calcu¬ 
late the radiation pressure (or scattering) force, F sca (r) = AQ sca wRI(r)/c, and absorption force, 
F abs (r) = 4Q^wRI(r)/c, exerted on cylindrical nanowires, where w is the beam waist, R the 
radius of the nanowire, and I{r) is the intensity distribution of the incident laser light along the 
width of the nanowire. Qabs is the absorption coefficient calculated using scattering and extinction 
efficiencies 0. These forces are plotted in Suppl. Fig. [TJb). It is evident that the absorption force 
is much weaker than the scattering forces. Experimentally, we have observed that the bolometric 
force F bo i is the dominant photo-induced force on the nanowire, as the Au/Ag 2 Ga bimorph mate¬ 
rial responds very efficiently to optically induced thermal bending. In our system, we estimate the 
bolometric force to be around 100 times bigger than the radiation pressure force. 

The main detection method involves the use of scattered light in an interferometric manner. 
Displacement of the nanowire along the optical axis modulates the phase of the scattered light. 
The reflection interferes with a local oscillator in a homodyne measurement setup, providing a 
way to detect the phase modulation. Using a PID controller, we lock the phase of the homodyne 
system to the phase quadrature of the scattered light. Outside the control bandwidth, the motion 
of the nanowire is then proportional to the error signal of the control loop. Light interacting with 
the nanowire and transmitted through the microscope lenses is also used for detection with a split 
detector. This method detects the motion of the resonator on the direction transverse to the optical 
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axis and is complementary to the interferometric measurement. A juxtaposition of the two methods 
is presented in Suppl. Fig. |TJc). Here we see two modes of vibration at different frequencies. In 
fact these peaks are due to spatially orthogonal modes of the nanowire. 

The homodyne measurement efficiency of a vibrational mode depends on the angle between 
the direction of oscillation and the axis of detection: if the oscillation is perfectly aligned with the 
optical axis, for example, that mode will be fully detected by the homodyne system and vibrations 
that are spatially orthogonal will not be observed. An example of this situation is shown in Suppl. 
Fig. [2]^ a) where we have rotated the nanowire until only one mode is visible. When the detection 
beam is incident with a 45° angle relative to the two orthogonal modes, both can be detected. This 
is the case in Suppl. Fig. (TJc) where both the split detector and the homodyne system can see the 
motion of the nanowire. The angle in this case is not exactly 45° since the peak heights are not 
identical for both methods, rather the largest peak for the homodyne signal is the smallest in the 
split detector signal, and vice-versa, as we expect. 

Supplementary Note 2 - Deflection. When subject to feedback from an intensity-modulated 
beam, the raise in bulk temperature of the nanowire, taking into account the reflectivity of the 
gold layer [4:] and an amplitude modulation depth of 0.1% for a 1 mW laser beam, is estimated 
to be around 10 K. It has been demonstrated that such temperature increase can cause thermally- 
induced nanomechanical deflection of a few nm of hybrid nanowires [5J. 

We were able to estimate the deflection by systematic modulation of the driving field for different 
powers. By turning the field on and off at low frequency to produce a variation in the detection 
signal, and by calibrating the amplitude of these variations to the total amplitude of the interference 
fringes of the unlocked homodyne signal (equivalent to one wavelength), we were able to measure 
the thermal time response (Suppl. Fig.[3]^a)-(b)) and observed deflections of about 60 nm with less 
than 5 mW of power (Suppl. Fig. [3](c)). These results reveal that a minimum power is required to 
drive the nanowire. 

The sign of the feedback gain can determine whether both modes can be cooled or heated simulta¬ 
neously. If the incident beam (or, equivalently, the nanowire) is rotated by 90°, the relative phase 


18 


between one of the modes and feedback will be opposite, while the other will be the same as be¬ 
fore. This will produce cooling of one mode and heating of the other, as shown in Suppl. Fig.[2}b). 
A rotation of a further 90° will again sync the phase of the two modes relative to the feedback, 
although now the phase that was originally used for cooling will produce heating, and vice versa. 
The full-rotation periodicity is a consequence of the thermal nature of the force inducing direc¬ 
tional bending, which depends purely on geometrical properties and not on the direction of the 
incident light: if that were to be the case, the behaviour after a 180° rotation should be the same as 
that of the nanowire in the original position, as the movement caused by radiation pressure forces 
is always in-phase with the driving field. 


Supplementary Note 3 - Feedback cooling. In general, the equation of motion [|6j of the nanowire 
is given by 


x + 7mT + uj^x = F lh + y / F n (t')h n (t - t')dt ', 


(3) 


where cu m is the mechanical frequency, 7 m is the mechanical damping, F lh is the thermal force 
driving the Brownian motion, and F n includes all of the other forces such as radiation pressure 
and photothermal forces. The term h n (t) = 1 — e~ t/Tc denotes the response time of the structure 
for different forces with time constant r c . In the case of the radiation pressure force this is al¬ 
most instantaneous (h lv (t) ~ 1), but response time of photothermal forces depends on the thermal 
conductivity of the material and geometry of the structure. The response time of nanocylinders 
to temperature variations depends on radius and thermal diffusivity k as r c ~ R 2 / (4k) Ma¬ 

terial properties and size strongly affect the thermal diffusivity, and for nano-scale objects it can 
be almost 2-3 orders of magnitude smaller compared to the bulk material [8]|9j due to phonon 
scattering overcoming the phonon-phonon coupling. Assuming an effective thermal diffusivity of 
approximately 10~ 6 m 2 s -1 , the thermal relaxation time of the nanowire can be on the order of 
ns or smaller, implying a response faster than its mechanical vibrations. The heat generated by 
the laser power will diffuse inside the nanowire more rapidly than the modulation, allowing the 
nanowire to react quickly to temperature fluctuations. 

To explain some of the observed effects, such as the relative behaviour between orthogonal modes, 
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and to be able to predict and have direct control on the system, we developed a simple model for 
our feedback setup that takes into account the transfer function of the system. The model includes 
a stability analysis of the feedback loop, from which it emerges that in our parameter regime we 
are far from any singularity. Thanks to the comparison of the total transfer function for both phase 
and gain with the experimental data collected we gained valuable insight on the spatial orientation 
of the modes. Suppl. Fig. [4] shows the experimental transfer functions and simulated feedback re¬ 
sponses for the nanowire at two different orientations. The experimental transfer function is used 
to infer the orientation angle of the nanowire. Using this information, the model predicts accurately 
the behaviour of the two orthogonal modes and their relative detection. 


Supplementary Note 4 - Noise squashing. For high gains, the measurement noise becomes sig¬ 
nificant and it is important to include a correction term in the expression for the displacement 


power spectrum observed [10|, 
2 7m T 0 


S.r.r ( -C ) 


m 2 (o;2 - cc 2 ) 2 + (1 + g) 2 ^u 2 


+ S, 


det" 


~ ^ 2 ) 2 + 7> 2 

ul-u 2 ) 2 + (1 + c/) 2 7m^ 2 ’ 
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which depends on the spectral density of detection noise, Sd et . The temperature inferred from this 
spectrum is 


rp _ T 0 g 2 mul y m 

-Off — Z-, -r —- — -Odet, 


1 + g 1 + g 4 k 


(5) 


B 


presenting now a lower bound T min = y '"y^' 7 ' 1 ,S det . 

A characteristic of these correction terms is that the displacement spectrum can get lower than the 
shot noise level. This is in fact typical of feedback systems, since the noise measured by our main 
detection scheme is inside the feedback loop and is not a faithful representation of the real noise of 
the system. The random fluctuations on the homodyne local oscillator beam are part of the infor¬ 
mation sent through the feedback to the nanowire, and are transferred to the phase of the homodyne 
detection beam once it is scattered back by the oscillator. When the scattered light interferes with 
the local oscillator, the noise on the two beams is now correlated and can destructively interfere. 
As a result, the measured noise can appear squashed below the shot noise level CD; out-of-loop 
detection would reveal this to be a property inherent to the feedback. Other experiments have 
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also observed this phenomenon [TOj. To deliberately observe noise squashing, we increased the 
local oscillator power to have larger shot noise. As it can be seen in Suppl. Fig. [5][a), applying 
the feedback with significant gain results in noise squashing where the amplitude noise measured 
by homodyne goes below the shot noise. The minimum temperature T min is reached when the 
Lorentzian profile lies flat on the noise level; for higher gain, squashing becomes manifest and 
the effective temperature gets higher again. An example of inferred temperature for increasing 
powers of the cooling beam (i.e. increasing g) is shown in Suppl. Fig. |5Jb), where the minimum 
temperature is reached at ~ 0.5 mW of power. 


Supplementary Note 5 - Estimation techniques. Given the dynamics of some system, the (ex¬ 


tended) Kalman filter [ 121 processes a series of noisy measurements of the system in order to keep 
track of a statistically optimal estimate of its underlying state. Our implementation of the extended 
Kalman filter was based on the code at reference [ 13 ]. Since the specifics of the filter are quite 
involved we will not discuss them in full detail here, and only indicate the general procedure. For 


a more in-depth discussion we refer the reader to the references [13] (for implementation details) 
and [ 14-fl6| (for complete introductory discussions). 


In our implementation we store the system state as a 4-vector containing the oscillator’s 
position x, velocity x, damping ratio 2( = x m /oj m and mechanical frequency u m . We also keep 
track of a covariance matrix which describes the uncertainty in this state vector. At each timestep 
(when a new measurement is taken, which is every 2 ps in our case), the filter updates in two stages. 
In the first “prediction” stage, the initial value problem describing the evolution of the system from 
the current state is solved using the Runge-Kutta method (RK4) to predict the state at the next 
timestep. In the RK4 approximation, the next value is determined by the present value plus the 
weighted average of four increments, where each increment is the product of the time interval 
dt = 2 ps with the derivative of the vector FT) = (x(t), x(t), £, c u m ) (calculated from the equation 
of motion for the system). 


In the second “updating” stage, the predicted oscillator position is compared to the measured 
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value, and the estimated state vector is updated to reduce this difference in accordance with the 


relative uncertainties (that is, more weight is given to either the measurement or the estimated state 
depending on which quantity has the lower uncertainty). The resulting state vector estimate is used 
as the initial state for the next timestep, and so on. This procedure yields an estimated trajectory 
of the oscillator, and from this we may extract a signal as described in the main text. 

To apply this algorithm we needed to define initial values for the mechanical frequency 
6o’ m /2 tt and linewidth 7 m /27r. These were first roughly estimated from the raw data, and then 
repeatedly adjusted in order to maximise the resulting SNR of the filtered trajectories. The final 
values we used were u m /2iT = 339.9 kHz and 7 m /27r = 0.53 kHz. We note that this linewidth 
is smaller than that estimated from the raw data (roughly 0.8 kHz). It is not entirely clear why 
it is advantageous (in terms of extracting the SNR) to model the oscillator with a longer thermal 
relaxation time, but a possible explanation is obtained by considering the behaviour of the Kalman 
filter’s prediction algorithm. In particular, we recall that to determine the SNR we turn off the 
updating stage of the filter after 1 ms, and quantify the signal by comparing the predicted trajectory 
to the observed trajectory. The predicted trajectory is simply a sinusoid decaying with the thermal 
relaxation time, whereas in reality these oscillations do not decay since the nanowire is always 
being driven by thermal noise (even when no signal is applied). A longer thermal relaxation time 
in the model causes the predicted oscillations to decay less slowly and thus, on average, provide a 
better prediction of the true oscillations. 

We implemented the virtual cooling method as described in reference GD- The idea is that 
the raw measurement record of the oscillator’s position (taken in the absence of feedback) may 
be post-processed to yield the measurement record that would have been obtained if a particular 
feedback scheme had been employed. We assume an ideal feedback force which is switched off at 
time to. 



( 6 ) 


where x is the measurement record, g is the feedback gain, H is the Heaviside step function and 5 
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is the Dirac delta function. As per |T7J, we then define a function 


/•OO 

h(t,r) := / 5(t' 


r)gH(t 0 - t')x(t - t')dt' = gH(t 0 - r)x(t - t), (7) 


J — OO 


where x is the mechanical susceptibility. The function h describes the response of the system at 
time t to feedback resulting from the position of the oscillator at time r, and is used to simulate 
the feedback process. Specifically, it is shown in reference G3 that the simulated measurement 
record x * satisfies the Fredholm equation 



( 8 ) 


To solve this equation we again follow the procedure of reference [17] and discretize to 1000 
timesteps of length 2gs, so that h is represented by a 1000 x 1000 matrix HI, while x and x* are 
represented by vectors x and x t of length 1000. Equation [8] becomes x* — Hx, = x, which is 
solved to obtain the simulated measurement record x* = (I — HI) _1 x (where I is the 1000 x 1000 
identity matrix). 

As with the Kalman filter, the virtual cooling method required several initial parameters to 
be defined. We determined these as follows. First, rough values of mechanical frequency u m /2ir, 
linewidth 7 m /27r and time t 0 were estimated from the raw data, and then the gain g was varied 
to maximise the SNR of the resulting simulated measurement record. The first three parameters 
were then repeatedly adjusted in order to maximise the peak SNR obtainable by varying the gain. 
The final values used were cc m /27r = 339.722 kHz, 7 m /27r = 0.85 kHz, t 0 = 0.896 ms and g = 
-1.6784 x exp(— 0.00004*) x 10~ 9 . 

These two estimation techniques allow a raw trace to be filtered in order to enhance the SNR 
of any signal present in the oscillations. We denote by // S nr this enhancement factor (or rather, 
the average of enhancement factors over all traces). Note that it does not make sense to discuss 
this factor in relation to feedback cooling, since if we are given a cooled trace we cannot “undo” 
the cooling to calculate the SNR that we would have obtained without cooling. Thus for feedback 
cooling we may define the enhancement in the average SNR (by comparing the average SNR with 
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feedback cooling to the average SNR without), but not the trace-by-trace enhancement factor ?/ SNR , 
therefore it is not reasonable to estimate the standard deviation of enhancement factor for feedback 
data. 
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Supplementary Figures 
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Supplementary Figure 1: Scattering simulation and detection, (a) Total nanowire scattering 
efficiency as a function of its radius for two different wavelengths, A = 780 nm and 1064 nm. Inset 
shows the total scattering efficiency in the XY-plane (seen from below) around the nanowire placed 
at the origin, when laser light is incident from right to left, (b) Absorption and scattering forces 
exerted on the nanowire as a function of its radius for A = 1064 nm and input power of 1 mW. 
(c) Detection of two orthogonal modes of the nanowire, using interferometry of the reflected beam 
(dark) and split detection of the transmitted beam (light). The two peaks are detected with different 
efficiency depending on the angle of the mode with the axis of detection, parallel to the optical axis 
for reflection and perpendicular to it for transmission. 
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Supplementary Figure 2: Effect of detection angle and phase on cooling, (a) The displacement 
spectrum of the nanowire used for cooling as shown in the main article, with (blue) and without 
(red) the feedback signal. The nanowire has been rotated so that only a single vibrational mode is 
visible in this frequency range. Note that this is 45° rotated compared to the situation for the data 
shown in Fig. 2 of the main manuscript. The solid lines represent the theory, (b) The displacement 
spectrum of the two modes of the same nanowire rotated by a further 45° relative to its orientation 
in (a). By swapping the sign of the gain, the behaviour of the feedback is altered in a mode selective 
manner. The red trace is obtained without feedback, blue and green data correspond to feedback 
with 0 and 7r phase shift, respectively. 
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Supplementary Figure 3: Characterization of the nanowire deflection, (a) Deflection of the 
nanowire measured using the homodyne signal while power is modulated at a rate of 17 kHz 
(shown as a square wave), (b) Exponential fit of deflection amplitude decaying after the driving 
force is switched off. The high-frequency oscillations on top of the exponential decay represent 
oscillations at the mechanical frequency, (c) The quasi-static displacement measured at different 
laser powers. The line is a guide to the eye. We infer a bolometric force of about 10 -10 N at 1 mW 
of power, taking into account a spring constant of 2 x ICE 2 N m _1 . 
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Supplementary Figure 4: Transfer function and feedback modelling for opposite orientations. 

(a) Experimental transfer function of gain (green) and phase (red) for the nanowire in the orien¬ 
tation used for the data shown in the main article, (b) Feedback model prediction for the transfer 
function at different gain values. Increasing gain causes heating of both modes, while gains at 
opposite phase lead to cooling as shown in the main text. The dashed line indicates absence of 
feedback, (c) Similar to (a), for the nanowire rotated by 90°. (d) Feedback model prediction 

for the new orientation inferred by the transfer function. The two modes are out of sync, and one 
is heated while the other one is cooled. For comparison, Suppl. Fig. |2|b) shows the experimental 
power spectrum corresponding to this orientation. 


29 



























(a) (b) 



Supplementary Figure 5: Noise squashing and feedback cooling at higher shot noise levels, (a) 
Noise spectrum showing the thermal amplitude noise of the nanowire (red), shot noise (black) and 
feedback result (blue) representing the noise squashing, (b) Temperature of the first vibrational 
mode of the nanowire as a function of power of the cooling beam. The solid lines in (a) and (b) are 
theory fits. 
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